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ABSTRACT 

Co-channel  and  adjacent-channel  interference  predictions  for  the 
Instrument  Landing  System  are  presented  in  this  report.     Information 
on  the  various  types  of  ground  facilities  involved  is  given  first.     Propa- 
gation mechanisms  applicable  to  the  108  to  112  Mc/s  frequency  range 
are  "discussed  second  together  with  the  calculation  of  transmission  loss 
and  its  variability.     Third,    the  statistical  treatment  of  interference 
problems  is  explained.     Finally;    the  results  of  the  study  are  presented 
as  curves  of  normalized  desired-to-undesired  signal  strength  ratios 
versus  distance  from  the  desired  station.     Aircraft  altitudes  of  1,  000, 
6,  250,   12,  000,   and  18,  000  feet  along  with  station  separations  ranging 
from  20  to  330  nautical  miles  were  considered. 

Detailed  procedures,   mathematical  formulas,   and  computer 
programs  used  are  discussed  in  the  Appendices. 

KEY  WORDS 

Instrument  Landing  System 
Interference  Predictions 
Propagation 


INTERFERENCE  PREDICTIONS  FOR  THE 
INSTRUMENT  LANDING  SYSTEM 


by 


G.    D.    Gierhart  and  M.    E.    Johnson 


1.     Introduction 

Increasing  air  traffic  density  together  with  fast,    high-flying  jets 
have  made  the  use  of  reliable  air  navigation  aids  more  important  than 
ever  before.     In  expanding  the  present  complex  of  navigation  aids  to 
meet  future  demands,    consideration  must  be  given  to  potential  inter- 
ference between  facilities  operating  on  the  same  or  on  adjacent  channels. 
The  amount  of  interference  is  a  function  of  the  desired-to-unde sired 
signal  ratio  at  the  aircraft  antenna  terminals;  as  both  signals  vary  with 
time  and  aircraft  location,    the  ratio  varies  as  well,    and  the  interference 
becomes  dependent  on  time  and  location.     Because  of  the  nature  of  radio 
wave  propagation  in  the  frequency  ranges  used,    the  variations  of  the 
received  signals  and  of  the  interference  ratios  are  best  described 
statistically.     The  large  number  of  possible  conditions  dictates  the  use 
of  a  digital  computer  with  programs  that  take  into  account  all  variables 
as  well  as  the  fixed  equipment  parameters. 

The  air  navigation  aid  treated  in  this  report  is  the  Instrument 
Landing  System  (ILS).     This  system  includes  a  runway  localizer,    a 
glide  path,    and  marker  beacons.     However,   in  this  study  only  the  local- 
izer is  considered,    since  it  is  most  susceptible  to  co-channel  and 
adjacent-channel  interference.     The  ILS  localizer  operates  in  the  108  to 
112  Mc/s  frequency  range  and  shares  this  range  with  the  VHF  Omnirange 


(VOR)  in  such  a  way  that  VOR  facilities  have  to  be  considered  as  the 
source  of  adjacent-channel  interference. 

At  very  high  frequencies  in  the  30-300  Mc/s  range  (VHF),    propa- 
gation of  radio  frequency  energy  is  affected  by  the  lower  atmosphere 
(the  troposphere),    specifically  by  variations  in  the  refractive  index  of 
the  atmosphere.     The  terrain  along  and  in  the  vicinity  of  the  great  circle 
path  between  transmitter  and  receiver  also  plays  an  important  part. 

Within  the  last  decade  a  number  of  methods  and  procedures  have 
been  developed  to  calculate  field  strength  and  its  variability  at  VHF. 
In  the  work  discussed  here  procedures  were  followed  which  have  been 
used  by  the  Central  Radio  Propagation  Laboratory  to  predict  in  a 
statistical  manner  the  effect  of  terrain  and  atmosphere  on  the  variability 
of  field  strength  and  on  the  performance  of  radio  systems  [Rice,    Longley, 
Norton,    and  Barsis,    196 5 J  .     It  is  also  convenient  to  use  the  concept  of 
transmission  loss  [Norton,    1953  and  1959]  ,    which  is  the  ratio  of  power 
radiated  to  the  power  that  would  be  available  at  the  receiving  antenna 
terminals  if  there  were  no  circuit  losses  other  than  those  associated 
with  the  radiation  resistance  of  the  receiving  antenna.     It  is  used  under 
the  assumption  that  losses  in  the  antenna  circuits  can  be  neglected,    and 
it  is  usually  expressed  in  decibels.     Methods  used  for  its  calculation  as 
a  function  of  path  length,    terminal  heights,    and  carrier  frequency  are 
discussed  in  Section  3.     Computation  techniques  are  discussed  in 
Appendix  B. 

After  some  initial  calculations,    parameters  for  the  ILS  localizer 
and  the  VOR  were  assembled  into  a.  computer  program  which  yields 
normalized  signal-to-interference  ratios  for  given  probability-of- 
service  values  as  a  function  of  aircraft  location  in  relation  to  the  desired 
and  the  unde sired  ground  stations.     The  normalization  involves  the 
assumptioii  that  both  the  power  and  antenna  gain  associated  with  the 
desired  and  unde  sired  station  are  the  same. 


Results  in  the  form  of  curves  showing  these  ratios  for  both  co- 
channel  and  adjacent-channel  cases  were  obtained.     Section  5  contains 
the  information  required  to  interpret  the  normalized  curves  in  terms  of 
parameters  associated  with  three  types  of  ILS  ground  equipment. 

2.     System  Parameters 

System  parameters  are  based  on  several  types  of  ILS  ground 
equipment,    a  standard  VOR  ground  installation,    and  passenger -jet 
aircraft  antennas.     Characteristics  of  three  ILS  localizers  are  listed 
in  Table  1.     Other  equipment  exists,    but  consideration  of  these  three 
is  sufficient  for  practical  purposes. 

TABLE  1 
Characteristics  of  ILS  Localizers 

Standard  Directional  Low  Cost 

Carrier  Power*  +  20  dBW  +  20  dBW  +  10  dBW 
Array  Type                                          8 -Loop                     V-Ring  V-Ring 

Antenna  Gain*  +  4  dB  +  12  dB  +12  dB 
Array  Height                                      5.  5  ft                         7.5  ft  7.  5  ft 

In  Table  1  carrier  power  refers  to  the  carrier  power  delivered  to 
the  terminals  of  the  carrier  antenna,    and  antenna  gain  refers  to  the 
main  lobe  free  space  gain  of  the  carrier  antenna  with  reference  to  an 
isotropic  radiator.     Figures  1  and  2  show  relative  gain  as  a  function  of 
azimuth  angle  for  the  carrier  portion  of  the  8 -Loop  [Civil  Aeronautics 
Administration,    1957]    and  V-Ring  arrays,    respectively. 

*See  text. 
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Figure   1.       8-Loop  Array  Antenna  Pattern 
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Figure    2         V -Ring  Array  Antenna  Pattern 


Both  the  ILS  localizer  and  the  VOR  use  horizontally  polarized 
antennas  and  were  assumed  to  be  operating  at  110  Mc/s.     Other  param- 
eters used  for  the  VOR  are  listed  in  Table  2. 

TABLE  2 
VOR  System  Parameters 

Carrier  Power  Delivered  to  Antenna  20  dBW 

Ground  Antenna  Type  4-Loop  Array  (located 

above  counterpoise) 

Maximum  Antenna  Gain  Relative  to 

Isotropic  Antenna  2. 15  dB 

Free  Space  Horizontal  Pattern  Approximately  Circular 

Free  Space  Vertical  Pattern  Approximately  Like  Dipole 

Counterpoise  Diameter  52  feet 

Antenna  Height  Above  Counterpoise  4  feet 

Counterpoise  Height  Above  Ground  12  feet 

Transmission  line  loss  associated  with  the  airborne  terminal  was 
considered  to  affect  both  the  desired  and  the  undesired  signals  equally 
and  was  neglected. 

Aircraft  antenna  gain  statistics  were  obtained  from  modeling 
study  data  based  on  an  E -cavity  type  VOR  antenna  in  the  vertical  stabi- 
lizer of  passenger  type  jet  aircraft  [Convair,    1959;  Commercial  Jetstar, 
1959]  .     Only  the  forward  ±  20  °  of  azimuth  were  considered  in  obtaining 
statistics  for  gain  toward  the  desired  station,    but  two  sets  of  statistics 
were  developed  for  gain  in  the  direction  of  the  undesired  station.     One 
set  considered  only  the  rear  ±  20°  of  azimuth  and  the  other  considered 
all  azimuth  angles  as  equally  likely.     From  these  statistics  a  single 
distribution  was  established  for  the  ratio  of  antenna  power  gain  in  the 
direction  of  desired  station  to  that  in  the  direction  of  the  undesired 


station.     This  ratio,    expressed  in  decibels,   is  denoted  by  the  symbol 

R        and  the  cumulative  distribution,     R     (p)  ,      used  to  account  for  air- 

craft  antenna  gain  is  shown  in  figure  3.     Also  shown  in  figure  3  are  two 

additional  cumulative  distributions  of  R  .      that  resulted  from  the  above 

A 

mentioned  analysis  and  were  used  as  a  guide  to  establish  the  R        distri- 
bution  used  in  the  calculations. 


3.     Transmission  Loss  Calculations 

Figure  4  shows  a  typical  configuration  of  an  aircraft  (representing 
the  receiving  terminal),    a  desired  navigational  transmitting  facility, 
and  an  unde sired  navigational  transmitting  facility.     All  three  are 
aligned  along  a  great  circle  path,    and  for  simplicity  assumed  to  be 
above  a  smooth  surface.     In  the  example  drawn,    the  aircraft  is  within 
the  radio  horizon  of  the  desired  facility,    but  beyond  the  radio  horizon  of 
the  interfering  station.     The  distances  along  the  great  circle  path  from 
a  point  vertically  below  the  aircraft  to  the  desired  and  the  unde  sired 
station  are  denoted  by  d       and  d      ,     respectively.     The  aircraft  is  at 
a  height    h       above  the  terrain.     The  angle     6    between  the  horizon  rays 
from  the  aircraft  and  the  interfering  station  is  an  important  parameter 
in  the  calculation  of  transmission  loss  for  beyond -the -horizon  paths 
[Norton,    Rice,    and  Vogler,    1955].     For  the  assumption  of  a  smooth 
spherical  earth,      9    is  also  the  distance  between  radio  horizons  expressed 
in  angular  measure,    and  it  is  therefore  called  the  "angular  distance,  " 
as  indicated  on  figure  4.     As  is  customary  in  the  analysis  of  tropo- 
spheric  propagation  paths,    first  order  allowance  for  the  refractive 
effects  of  the  atmosphere  was  made  by  assuming  an  earth  radius  4/3 
times  the  actual  radius.     This  permits  radio  rays  to  be  drawn  as 
straight  lines,    and  simplifies  the  determination  of  geometric  parameters. 

Transmission  loss  calculations  were  accomplished  by  (a)  calcu- 
lating a  reference  value  of  basic  transmission  loss,    (b)  calculating  a 
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cumulative  distribution  for  long-term  variations,    (c)  calculating  a 
cumulative  distribution  for  short-term  variations,    and  (d)  calculating 
the  cumulative  distribution  of  transmission  loss  by  combining  the 
results  of  previous  calculations.     More  detail  on  these  four  steps 
follows: 

(a)     A  reference  value  of  basic  transmission  loss  was  calculated 
in  accordance  with  methods  given  by  Rice,    et  al.    [1965]  . 
This  reference  value,    as  explained  by  Norton,    et  al.    [1955]  , 
is  the  median  of  all  hourly  medians  within  "Time  Block  2,  " 
representing  the  hours  1:00  p.m.    to  6:00  p.m.,    November 
through  April.     In  calculating  "basic  transmission  loss"  the 
antennas  are  considered  to  be  isotropic  with  unity  power  or 
voltage  gain  in  all  directions.     Actual  antenna  characteristics 
are  brought  in  later  on  when  converting  from  basic  trans- 
mission loss  to  transmission  loss.     Within  the  radio  horizon, 
basic  transmission  loss  was  calculated  using  geometric  optics 
methods,    including  interference  between  the  direct  and  the 
ground-reflected  ray.     For  the  desired  ILS  localizer  and  VOR 
propagation  models,    specular  reflection  was  assumed  with 
reflection  coefficients  equal  to   -0.9.      Because  of  the  irregular 
nature  of  the  terrain  (buildings,    etc.    included)   surrounding 
the  unde  sired  ILS  stations,    specular  reflection  was  considered 
less  dominating,    and  a  combination  of  specular  and  diffused 
reflection    was    assumed.      The  primary  effect  of  this  assump- 
tion was  to  lower  transmission  loss  values  associated  with  the 
undesired  ILS  relative  to  comparable  values  for  the  desired 
ILS.      This  is  discussed  further  in  Appendix  A.      Beyond  the 
radio  horizon,    basic  transmission  loss  was  calculated  using 
smooth-earth  diffraction  or  forward  scatter  models,    depending 


10 


on  the  path  distance  involved.     Calculations  for  both  models, 
and  the  method  of  properly  combining  diffraction  and  scatter 
fields  if  they  are  of  comparable  magnitude,    were  based  on 
procedures  contained  in  a  Technical  Note  by  Rice,    et  al.    [1965]. 

(b)  Long-term  variations  in  basic  transmission  loss  were  estimated 
for  a  continental  temperate  climate  by  means  of  the  time  varia- 
bility function    V(p,d).     This  empirical  function  was  calculated 
in  accordance  with  methods  given  by  Rice,    et  al.    [1965]  .     It 
was  used  to  determine  the  cumulative  distribution  with  time  of 
hourly  median  basic  transmission  loss  values  relative  to  the 
reference  values  calculated  under  (a)  above  for  all  hourly 
median  basic  transmission  loss  values  during  the  year  as 
functions  of  path  length,    angular  distance,    terminal  height,    and 
carrier  frequency. 

(c)  In  addition  to  the  distributions  of  hourly  medians  representing 
long-term  variations,    short-term  (usually  within -the -hour) 
distributions  of  the  received  signal  levels  had  to  be  estimated. 
Short-term  variations  in  this  particular  application  are  princi- 
pally due  to  two  causes.     One  is  the  inherent  short-term  fluctu- 
ation of  the  tropospheric  signal  ascribed  to  phase  interference 

of  rays  reflected  from  small  layers  or  scattered  from  refractive 
index  discontinuities,    or  to  reflections  from  ground  irregular- 
ities.    The  other  is  the  pattern  of  the  aircraft  antenna:     numerous 
small  lobes  cause  gain  changes  with  varying  bearings  which  can 
be  represented  by  a  cumulative  distribution  of  antenna  gain  with 
time  as  the  aircraft  moves  through  space.     However,    it  was 
more  efficient  to  neglect  the  effect  of  the  aircraft  antenna  gain 
in  the  initial  transmission  loss  calculations  and  include  it  in 
later  calculations  as  the    R.(p)     distribution  discussed  in 
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Section  2.     Short-term  fading  was  described  by  a  cumulative 
distribution  function    V     (p,  0)     which  is  based  on  fading  range 
data  given  by  Janes  [1955]  . 

(d)  To  obtain  the  cumulative  distribution  of  transmission  loss  the 
functions  discussed  in  a,  b,  and  c  above  were  combined.  The 
mechanics  of  combining  cumulative  distributions  is  discussed 
in  more  detail  in  Appendix  B. 

The  prediction  of  interference  conditions  requires  a  knowledge  of 
the  time  distribution  of  transmission  loss  or  field  strength  at  many 
points  in  space;  therefore,    time  variability  calculations  have  to  be 
performed  for  these  points  with  distance  from  the  desired  and  unde sired 
station,    aircraft  height,    spacing  of  desired  and  undesired  ground  facil- 
ities,   carrier  frequency,    and  antenna  patterns  as  parameters. 

It  should  be  noted  here  that  the  one -hour  period  taken  as  the 
dividing  line  between  long-term  and  short-term  variations  is  somewhat 
arbitrary.     It  is  convenient  in  view  of  the  available  empirical  time 
variability  functions. 

4.     Interference  Between  Two  Stations 

As  shown  by  figure  4,    both  the  desired  and  the  undesired  signals 
arrive  at  the  aircraft  over  propagation  paths  characterized  by  the 
distances    d        and    d      ,    respectively,    and  by  the  aircraft  height.     The 
distances  are  measured  along  the  great  circle  path.     Both  signals  vary 
with  time,    and  the  distributions  of  signal  levels  were  calculated  in 
accordance  with  the  procedures  outlined  in  the  preceding  section.     Then 
the  ratio  of  desired  to  undesired  signal  exceeded  for  given  percentage - 
of-time  values  at  a  particular  aircraft  location  was  determined. 

The  ratio  of  desired-to-undesired  signal  can  be  expressed  as  the 
decibel  difference  of  desired  and  undesired  signal  levels  and  was 
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obtained  from  calculated  transmission  loss  values  and  other  system 
parameters.     The  distribution  of  this  ratio  will  be  denoted    D/U(p), 
where  p  is  that  percentage  of  time  during  which  a  given  value  of  D/U 
is  exceeded.      By  virtue  of  the  aircraft  being  in  motion,    time  variations 
also  include  variations  in  space.     As  the  actual  time  distribution  of 
D/U    may  vary  from  installation  to  installation  because  of  terrain 
characteristics  and  other  factors  not  taken  into  account  in  this  analysis, 
the  time  availability    p    may  be  interpreted  as  an  expression  of  relia-. 
bility  for  a  typical  installation.     The  concepts  of  "prediction  uncertain- 
ty" and  "service  probability"  in  the  sense  defined  by  Barsis,    et  al. 
[1962]    were  not  used  explicitly.     It  should,    however,    be  understood 
that  there  is  an  uncertainty  associated  with  the  D/U  predictions  of  this 
study  and  that  this  uncertainty  will  increase  under  conditions  where  the 
assumed  propagation  models  become  less  valid. 

As  an  example,     D/U(95)   =  10  dB    means  that  for  a  typical 
installation  the  ratio  of  the  desired-to-undesired  signal  is  equal  to  or 
greater  than  10  dB  during  95%  of  the  time.     Values  of    D/U(95)     are 
associated  with  the  variables  used  in  the  calculations.     These  include: 
(1)  interference  type  (co-  or  adjacent-channel),    (2)  aircraft  altitude, 
(3)  station  separation,    and  (4)  aircraft  distance  from  the  desired 
station. 

To  obtain  the  time  availability  of  the  desired-to-undesired  ratio 
at  any  point  in  space  it  was  necessary  to  properly  combine  the  cumu- 
lative distributions  of  (a)  transmission  loss  from  the  desired  station, 
(b)  transmission  loss  from  the  undesired  station,    and  (c)  the  antenna 
power  gain  ratios,     RA(p).     Actually  these  calculations  resulted  in  a 
cumulative  distribution  of  a  normalized    D/U.     The  process  for  con- 
verting normalized    D/U    to  actual    D/U    values  is  discussed  in 
Section  5. 
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5.     Results  of  the  Study 

The  results  of  this  study  are  in  the  form  of  normalized  prediction 
curves.      These  curves  may  be  used  to  estimate  the  service  limitations 
imposed  on  ILS  installations  by  co-channel  and  adjacent-channel  inter- 
ference.    Other  limitations  to  service  such  as  thermal  or  man-made 
noise  at  the  receiver  and  self  interference  caused  by  reflections  from 
airport  structures  or  other  aircraft  were  not  considered  in  this  study. 
In  particular,    an  "acceptable"  desired-to-undesired  signal  ratio  does 
not  imply  that  the  desired  signal  is   strong  enough  for  operational  use. 

In  this  section  and  in  Appendix  B  values  of  normalized    D/U(95) 
will  be  denoted  by  the  symbol    N{D/U(95)}.     These  normalized  values 
were  calculated  for  the  condition  when  the  two  ground  stations  and  the 
aircraft  were  aligned  over  the  same  great  circle  arc.     This  was  men- 
tioned previously  in  Section  3  and  is  implied  by  the  inset  in  figures   5 
through  10.     However,     N{D/U(95)}    values  can  also  be  obtained  for 
cases  v/here  the  assumed  alignment  is  not  valid  by  properly  interpreting 
the  station  separation  shown  on  the  curves.     Regardless  of  the  shortest 
distance  between  the  ground  stations,    the  station  separation,     S,    shown 
on  the  curves  should  always  be  regarded  as  the  algebraic  sum  of  the 
distance  from  the  aircraft  to  the  desired  station,     d    ,    and  the  distance 
from  the  aircraft  to  the  undesired  station,     d     ;  i.e.,     S  =  d      +  d     . 

5.1.     Service  Limitations  Due  to  Co -Channel  Interference 

N{D/U(95)}    curves  for  aircraft  altitudes  of  6,250,    12,000,    and 
18,000  feet  are  shown  in  figures   5,    6,    and  7,    respectively.     Desired 
values  of    D/U(95)     may  be  converted  to  values  of    N{D/U(95)}    which 
can  be  read  from  the  curves  by  the  following  procedure. 

(a)     Determine  the  value  of  co -channel  station  combination  factor, 
C,    from  Table  3,    on  page  18. 
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FREQUENCY     110  Mc/s 
ALTITUDE        6,250   FEET 
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Figure  5.      Co-Channel  IL.S  Interference;   6,  250  feet 
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FREQUENCY     110  Mc/S 
ALTITUDE      12,000  FEET 


STATION    SEPARATION  ,S,  AS  LABELED 
95%    RELIABILITY 
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Figure  6.      Co-Channel  ILS  Interference;  12,  000  feet 
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FREQUENCY     110  Mc/s 
ALTITUDE      18,000  FEET 


STATION    SEPARATION  ,S,  AS  LABELED 
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Figure  7.      Co-Channel  ILS  Interference;  18,  000  feet 
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(b)  Determine  the  azimuth  angle,      a,    between  main  lobe  maximum 
of  localizer  carrier  antenna  at  the  unde sired  station  and  the 
aircraft. 

(c)  Using     a    and  the  antenna  pattern  (figure  1  or  2)  appropriate  for 
the  unde  sired  station  type,    determine  the  gain  factor,    G,    of  the 
unde  sired  localizer  carrier  antenna  in  the  direction  of  the 
aircraft. 

(d)  Calculate    N  {D/U(95)  }  using 


N{D/U(95)}  =  D/U(95)  -  C  +  G 


(1) 


TABLE  3 
Co -Channel  Station  Combination  Factor,    C 


Desired  Station  Type 


Undesired  Station  Type 


Standard 


Directional 


Low  Cost 


Standard 
Directional 
Low  Cost 


0  dB 

10.  5  dB 

0.  5  dB 


-10.  5  dB 

0  dB 

-10  dB 


■0.  5  dB 

+  10  dB 

0  dB 


The  values  of    C    in  Table  3  were  calculated  from  information 
given  in  Table  1  by  using  the  following  equation: 


C=PD-PU+AD-AU+HD-HU 


(2) 


where        P       =  desired  station  carrier  power  in  dBW, 

P      =  undesired  station  carrier  power  in  dBW, 
A       =  free   space  antenna  gain  referred  to  an  isotropic 
radiator  for  the  main  lobe  of  the  desired  station 
localizer  carrier  antenna, 
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A       =  antenna  gain  similar  to    A     ,    but  for  unde sired  station, 
H      =  height  gain  factor  for  desired  station, 

H      =  0  dB  for  array  height  of  5.  5  feet  (8-Loop). 

H      =  2.  5  dB  for  array  height  of  7.  5  feet  (V-Ring). 
H      =  height  gain  factor  similar  to  H    ,    but  for  unde  sired 

station    (H      =  0  dB  for  adjacent-channel  calculations). 

Values  of    G    read  from  figures  1  or   2  represent  the  gain  of  the 
unde  sired  localizer  carrier  antenna  in  the  direction  of  the  aircraft 
relative  to  the  gain  of  the  same  antenna  in  the  direction  of  the  main 
lobe  maximum.     Because  of  this,    values  of    G    are  always  negative 
and  D/U(95)     will  have  its  worst  (lowest)  values  for  a  particular  pattern 
when    G    is  0. 


5.  2.     Service  Limitations  Due  to  Adjacent-Channel  Interference 

Adjacent-channel    N{D/U(95)}    curves  for  aircraft  altitudes  of 
1,  000,    6,  250,    and  12,  000  feet  are  shown  in  figures  8,    9,    and  10,    respec- 
tively.    Values  of    D/U(95)     maybe  converted  to  values  of    N{D/U(95)} 
which  can  be  read  from  the  curves  by  the  following  procedure: 

(a)  Determine  the  value  of  the  adjacent-channel  station  combination 
factor,     C,    from  Table  4. 

(b)  Calculate    N{D/U(95)}    using 

N{D/U(95)}  -  D/U(95)   -C  (3) 

TABLE  4 
Adjacent-Channel  Station  Combination  Factor,    C 

Desired  Station  Type 

Standard  Directional  Low  Cost 

1.85  dB  12.  35  dB  2.  35  dB 
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FREQUENCY     110  Mc/S 
ALTITUDE         1,000   FEET 


STATION    SEPARATION  ,S,  AS  LABELED 
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Figure  8.     Adjacent-Channel  IL.S  Interference;  1,  000  feet 
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FREQUENCY     110  Mc/S 
ALTITUDE       6,250   FEET 


STATION    SEPARATION  ,S,  AS  LABELED 
95%    RELIABILITY 
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Figure  9.      Adjacent-Channel  ILS  Interference;   6,  250  feet 
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FREQUENCY     110  Mc/S 
ALTITUDE      12,000   FEET 


STATION    SEPARATION  ,S,  AS  LABELED 
95%    RELIABILITY 
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Figure  10.     Adjacent-Channel  LLS  Interference;  12,  000  feet 
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The  values  of    C    in  Table  4  were  calculated  from  information 
given  in  Tables  1  and  2  by  using  the  equation  given  in  Section  5. 1. 
(H      =0  dB    for  VOR.  ) 

In  the  adjacent-channel  case  the  calculation  of    D/U(95)    is  some- 
what simpler  than  in  the  co -channel  case.     This  is  because  the  adjacent- 
channel  unde sired  station  is  always  a  VOR  and  the  VOR  has  an  omni- 
directional radiation  pattern  in  the  horizontal  plane;  thus  no  gain  factor 
is  required. 
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APPENDIX  A 
Propagation  Models 

A  smooth,    spherical  earth  model  was  used  to  calculate  trans- 
mission loss  for  the  ILS  and  VOR  system  at  110  Mc/s.     A  linear - 
gradient  atmosphere  was  assumed  for  the  initial  calculations,    so  that 
for  a  first  approximation,    radio  rays  could  be  considered  to  be 
straight  lines  above  an  earth  having  an  effective  radius    k    times  its 
actual  value.     In  conformance  with  usual  practice,    the  effective  earth 
radius  factor,     k,    was  assumed  to  be  4/3  for  the   "standard"  atmosphere. 
Calculation  methods  are  based  on  material  contained  in  a  Technical 
Note  by  Rice,    et  al.    [1965]. 

It  was  found  convenient  to  include  only  the  gain  of  the  transmitting 
antenna  in  the  definition  of  transmission  loss  as  used  in  these  calcula- 
tions.    Usually,    the  term  "transmission  loss"  includes  the  antenna 
gains  at  both  path  terminals,    whereas  the  term  "basic  transmission 
loss"  excludes  all  antenna  gains  [Norton,    1953  and  1959]  .     In  this  study, 

an  hourly  median  transmission  loss,     L      ,    is  defined  to  include  the  free- 

m 

space  gain  of  the   transmitting    antenna  as  a  function  of  the  elevation 

angle  as  well  as  the  effect  of  reflections  from  the  counterpoise  (VOR 

only)  and  from  the  ground. 

Within  the  radio  horizon,    values  of  the  hourly  median  transmission 

loss,      L      ,    were  calculated  using  geometrical  optics  methods  [  Kirby, 
m 

Herbstreit,    and  Norton,    1952]  .     These  methods  take  into  account  the 

interference  between  the  direct  and  the  ground-reflected  ray.     Figure  11 

shows  the  geometry  for  this  ray  interference  problem,    and  defines 

many  of  the  symbols  used  in  the  analysis.     Using  this  geometry,    curves 

of    L        versus  distance  were  calculated  for  each  assumed  aircraft 
m 

altitude.     The  following  constants  and  relationships  were  used,    with 


24 


distances  and  angles  measured  in  statute  miles  and  radians,    respec- 
tively: 

a    .  5280  statute  miles  =  4/3  times  the  actual  radius  of 

the  earth 
h    .      .      .      .     height  of  the  transmitting  antenna  above  the  reflecting 

surface 
h         ...     aircraft  altitude  above  the  reflecting  surface 
|R  I    .      .      .     magnitude  of  the  reflection  coefficient 
c    .  phase  angle  of  the  reflection  coefficient  relative  to     TT, 

and  assumed  to  be  zero 
D  .  divergence  factor 

d    .  distance  along  the  great  circle  path  to  a  point  on  the 

reflecting  surface  below  the  aircraft 

r    .  length  of  the  direct  ray  path 

o 

r    +  r  .     length  of  the  reflected  ray  path 

d    .      .      .      .     distance  from  transmitting  antenna  to  reflection  point 

d  ...     distance  from  aircraft  to  reflection  point 

La 

i|»   .  grazing  angle  at  reflection  point 

The  first  step  was  the  determination  of  the  grazing  angle     ijj.     For  a 
given  value  of    d    the  following  equations  were  used: 

d  =  d1  +  d2  (4) 

V  =hi  ~di2/2a  (5) 

V    -h2"d22/2a  {6) 

taniJJ=h1,/d1  =h2'/d2.  (?) 

Equations  (5),    (6),    and  (7)  are  approximate  relations.     In  order  to 
utilize  the  electronic  computer  in  solving  the  above  expressions  for 
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Figure   11.     Geometry  for  Ray  Interference  within  the  Radio  Horizon 
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grazing  angle,      ijj,    an  iteration  method  was  used,    assuming  for  each    d 

different  values  of    d      until  the  equations  were  satisfied  within  the 

desired  limits  of  accuracy,    which  was     ±10         miles.     Final  values  for 

\\i  for  each  assumed  total  distance    d    were  then  used  in  calculating 

various  auxiliary  angles  and  quantities  to  be  used  in  ray-tracing  and  in 

the  final  calculation  of  the  transmission  loss     L. 

m 

Using  the  previously  calculated  value  of    ijj   other  geometrical 
parameters  shown  in  figure  11  were  calculated  as  follows: 

sin  (j>      =    (a  cos  ijj)/(a  +  h  )  (8) 


sin  cj>       =    (a  cos  4»)/(a  +  h   )  (9) 


>x    =    tt/2  -  iji-  ^  (10) 


02    =    n/2  -  i|i-  ct>2  (11) 


r2     =    (a  sin  Qj/  (sin  <J>2)  (12) 


90     =    91  +  92  W 


■0     -     /  (hx  -  h2)2  +  4(a  +  h1)(a  +  h2)[sin  (0Q/2)]  2         (14) 


dn  (3    =    [r2  sin  (  tt  -  2i|>)]  /rQ  (15) 


Y01    =    P-  +-0  (16) 


Vrl   =    -(+■+  Bx)  (17) 
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The  divergence  factor    D    was  calculated  using 


D     = 


(18) 


1    + 


2dld2 

a  d  tan  ip 


Voltage  gains    g        and  g        of  the  ground  antenna  relative  to  that 
of  an  isotropic  radiator  were  calculated  by: 


gQ1    =    1.  28  cos  vQ1 


(19) 


and 


g         -    1.  28  cos  v 
rl  rl 


(20) 


The  phase  angle     A,   in  radians,    between  the  direct  and  reflected 
rays  due  only  to  the  difference  in  ray  length  was  calculated  using 


2.42  f         h'  h   ' 
A     -  Mc     1       2 


(21) 


where    f , ,       =    110  Mc/s. 

Mc 

Then    L        was  calculated  from  the  following  formula: 
m 


L        =    36.  58  +  20  log.,  f   . 
m  °10    Mc 


(22) 


+  20  1og10r0-10  1og1()|g01     + 


D  |R|  grl]  2  -  2g01sriD  lRl  cos  A} 


In  the  above  equations  distances  and  heights  are  measured  in 

statute  miles  and  angles  measured  in  radians.     The  distance    d    was 

converted  to  nautical  miles  after    L        was  calculated. 
m 
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For  large  distances  and  low  angles  the  linear  gradient  atmosphere 
predicts  too  much  bending  [Bean  and  Thayer,    1959]  .     To  correct  for 
this,    a  ray  leaving  the  ground  station  antenna  at  the  same  angle,      y      , 
as  the  direct  ray  (calculated  using  a  linear  gradient  atmosphere)  was 
traced  through  an  exponential  atmosphere  until  it  reached  the  aircraft 
altitude.     The  great  circle  distance  below  the  ray  was  then  used  in  the 
final  calculations  instead  of  the  distance    d    first  assumed.     The  expo- 
nential atmosphere  used  corresponds  to  the  reference  atmosphere 
given  by  Rice,    et  al.    [1965]    with  a  surface  refractivity  of    N     =  301. 

Calculations  of  transmission  loss  beyond  the  radio  horizon  were 
made  by  combining  smooth-earth  diffraction  with  tropospheric  scatter. 
The  short-cut  method  of  Vogler  [1961]    was  used  for  diffraction  comput- 
ations.    Scatter  computations  and  methods  of  combining  the  two  mecha- 
nisms follow  the  work  by  Rice,    et  al.    [1965]  . 

Calculations  for  the  beyond-the -horizon  case  are  based  on  the 
same  exponential  atmosphere  used  for  within-the -horizon  computations. 
The  electrical  constants  of  the  ground  assumed  for  calculations  of  the 
diffracted  wave  were  a  dielectric  constant,      c    =  15,    and  a  conductivity 
value,      cr  =  5  millimhos  per  meter.     At  the  carrier  frequency  used, 
however,    the  effect  of  ground  constants  is  small. 

Finally,    from  a  large  number  of  computations  both  within  and 
beyond  the  radio  horizon,    a  set   of  curves  was  drawn  for  discrete  air- 
craft altitudes  showing  the  expected  median  transmission  loss  as  a 
function  of  the  great  circle  distance,     d,   (assuming  the  antenna  gain  of 
the  ground  station  to  correspond  to  a  half -wave  dipole).     The  construc- 
tion of  these  curves  involved  a  certain  amount  of  blending  in  the  vicinity 
of  the  radio  horizon  [Kerr,    1964]  . 

Long-term  variability,    short-term  fading,    and  the  effective  gain 
variations  of  the  airborne  antenna  were  treated  statistically.     Long- 
term  variability  was  described  by  the  time  availability  function  V(p,d). 
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Equations  and  curves  given  by  Rice,    et  al.    [1965]   were  used  to  calcu- 
late the  V(p,  d)    applicable  to  all  hours  of  the  year,    for  a  continental 
temperate  climate. 

Short-term  fading  was  described  by  a  cumulative  distribution 

function    V    (p,  9)     which  includes  the  summing  of  a  constant  vector  and 
F 

a  Rayleigh-distributed  vector.     The  angular  distance,     9,   is  defined  to 

be  the  angle  between  horizon  rays  from  the  transmitter  and  receiver 

[Norton,    Rice,    and  Vogler,    1955]  ,    and  is  defined  on  figure  4  of  the 

main  body  of  this  report.     An  auxiliary  distribution  function    V     (p,  9) 

F 

was  developed  using,  1)  fading  range  data  for  199  Mc/s  from  a  paper  by 
Janes  [1955]  to  express  fading  range  as  a  function  of  9,  and  2)  figure  6 
and  Table  1  from  Norton,  Vogler,  Mansfield,  and  Short  [1955]  to  spec- 
ify a  cumulative  distribution,      V     (p,  9),    for  the  fading  range  given  by 

F 

a  particular  value  of    9.     The  distribution  functions    V    (p,  9)     and 

F 

V     (p,  9)    are  related  by 
F 

VF(p,  9)     =    V'F(p,  9)  -  V'F(50,9).  (23) 

It  was  unnecessary  to  consider  a  different  aircraft  antenna  gain 
for  the  direct  and  reflected  rays  because  their  arrival  angles  are 
nearly  equal  as  would  be  more  apparent  if  figure  11  was  drawn  to  scale. 
This  made  calculations  considerably  simpler  since  the  effect  of  aircraft 
antenna  gain  on    N{D/U(95)}     could  be  considered  in  transmission  loss 
calculations  without  including  the  above  mentioned  gain.     The  method 
involved  is  outlined  in  Appendix  B  and  uses  the  distribution    RA(p) 
shown  in  figure   3. 

Distributions  of  transmission  loss  values  were  calculated  using 
the  following  equation: 

L(p,d)      =    L       -  [V(p,d)*V_(p,  9)]    .  (24) 

m  x 
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In  this  equation  all  terms  except    L         are  cumulative  distributions,    so 

m 

that  the  calculation  of    L(p,  d)     required  the  use  of  a  convolution  process 

which  is  indicated  by  the  operational  symbol    *  ,    and  will  be  further 

discussed  in  Appendix  B.     As  previously  implied,    in  calculating    L(p,  d) 

an  isotropic  aircraft  antenna  was  assumed.     When  the  fading  range  was 

zero  the  effect  of    V    (p,  9)     could  be  neglected  and    V     (p,  9)     was 

F  F 

dropped  from  equation  (24). 

In  general  the  above  discussion  applies  to  all  propagation  models 
used  in  this  study.     However,    more  information  is  required  to  make  the 
description  of  each  model  reasonably  complete.     Information  relevant 
to  particular  models  is  contained  in  the  remaining  parts  of  this  appendix. 

A.l    Desired  Station  ILS  Model 

The  propagation  model  used  for  the  desired  ILS  localizer  carrier 
system  is  described  accurately,    but  not  completely  in  the  first  part  of 
this  appendix. 

In  the  calculations,    the  ground  antenna  was  assumed  to  be  a  single 

Alford  loop  5.  5  feet  above  ground.     The  total  gain  of  the  carrier  antenna 

and  the  carrier  power  were  not  considered  in  the  calculation  of 

N{D/U(95)  }.     These  can  be  taken  into  account  by  the  procedure  given  in 

Section  6  for  converting    N{D/U(95)}    to    D/U(95). 

For  this  model  the  fading  range  values  on  which    V    (p,  9)     was 

F 

based  were  approximated  by 

0  for  9  <  0 
1  +  12.49  dB  for  9  <9  <1° 
13.4  dB  for  9  >  1°     . 
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A.  2    Undesired  Station  ILS  Model 

The  propagation  model  used  for  the  unde  sired  ILS  localizer 

carrier  system  was  the  same  as  that  used  for  the  desired  ILS  station 

with  respect  to  assumptions  made  about  the  ground  antenna  and  the 

methods  used  to  calculate    L        beyond  the  radio  horizon.     It  differed 

m 

in  the  method  used  to  calculate    L        within  the  radio  horizon  and  the 

m 

approximation  used  to  describe  fading  range. 

The  within-the -horizon    L         curves  for  the  unde  sired  ILS  were 

m 

developed  by  utilizing  the     L        curve  used  for  the  desired  ILS  along 

m 

with  a  similar  curve  developed  by  essentially  neglecting  the  ground 

reflections.     Both  curves  were  plotted  on  the  same  sheet  with  loss 

measured  along  the  ordinate  by  a  scale  linear  in  decibels  and  distance 

from  the  station  measured  along  the  abscissa  by  a  scale  linear  in 

nautical  miles.     The  unde  sired  ILS    L        curve  was  then  obtained  by 

m 

interpolating  a  smooth  curve  between  the  curves  previously  drawn  such 

that  the  new  curve  merged  with  the  desired  ILS    L        curve  in  the  vicin- 

m 

ity  of  the  radio  horizon  and  merged  with  the  "free  space"  loss  curve  at 
about  10  nautical  miles  from  the  station.     This  procedure  was  used  in 
order  to  estimate  the  effect  of  terrain  roughness  around  the  undesired 
ILS  and  to  avoid  making  quantitative  assumptions  about  the  effect  of 
terrain  roughness  on  the  specular  reflection  coefficient.     It  assumes 
that  the  specular  reflection  coefficient  is   -0.  9  at  maximum  within-the  - 
horizon  range  (low  grazing  angle)  and  that  as  the  range  is  decreased 
with  fixed  altitude  (higher  grazing  angle)  the  effect  of  terrain  roughness 
serves  to  reduce  the  absolute  value  of  the  specular  reflection  coefficient. 
A  random  reflection  coefficient  of  0.  3  [McGavin  and  Maloney,    1959] 
was  assumed  to  be  associated  with  the  rough  terrain  surrounding  an 
undesired  ILS.     The  effect  of  this  reflection  was  accounted  for  by  the 
short-term  fading  function    V     (p,9).     For  this  model  the  fading  range 
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values  on  which    V    (p,  9)     was  based,    were 

F 

4.  7  dB  for  9  <  9° 

1  +  12.  49  dB  for  0  <    9  <  1° 

13.  4  dB  for  9  >  1°. 

A.  3    Unde sired  Station  VOR  Model 

The  propagation  model  used  for  the  undesired  VOR  localizer  fol- 
lows the  description  given  at  the  beginning  of  this  appendix  except  in 
the  region  where  reflection  from  the  counterpoise  was  considered. 
Remarks  made  concerning  antenna  assumptions  in  the  discussion  of 
the  desired  ILS  model  are  applicable  to  the  VOR  model  except  for  the 
specific  value  of  antenna  height.     The  fading  range  values  used  in  the 

desired  ILS  model  to  determine    V     (p,  9)     were  also  used  in  the  VOR 

F 

model. 

Two    L     -versus -distance  curves  were  calculated  for  each  air- 
m 

craft  altitude  within  the  horizon.     One  curve  was  used  for  reflections 
from  the  counterpoise,    where  applicable,    and  the  other  for  reflections 
from  the  ground.     Both  curves  were  blended  in  the  region  where  the 
reflection  point  was  close  to  the  edge  of  the  counterpoise.     In  order  to 
facilitate  this  blending,    the  curves  were  calculated  so  that  they  extended 
beyond  their  region  of  validity. 

For  reflections  from  the  counterpoise,    calculated  for  elevation 
angles  greater  than  5  degrees,    the  effective  earth  radius  was  increased 
by  the  height  of  the  counterpoise  above  ground  so  that  a  sphere  through 
the  level  of  the  counterpoise  became  the  reflecting  surface  with  a  re- 
flection coefficient  of  -9.9  and  the  antenna  was  assumed  to  be  an  Alford 
loop  4  feet  above  this  surface.     The  aircraft  altitude  was  also  adjusted 
to  compensate  for  the  counterpoise  height.     In  these  calculations    d 
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was  very  much  smaller  than    d       and  the  following  two  approximations 
were  used  to  simplify  the  calculations: 


Pj    S   h^sin  i|i  (25) 

)      =    (r    cos  ljj )/a  .  (Z6) 


Beyond  the  region  where  reflection  from  the  counterpoise  was 
important    L         was  calculated  using  the  methods  outlined  at  the  begin- 
ning of  this  section.     The  antenna  was  assumed  to  be  an  Alford  loop 
located  16  feet  above  the  ground. 
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APPENDIX  B 
Computation  Techniques 

Calculations  necessary  to  produce  all  desired  distributions  of 
transmission  loss  values  or  power  levels  as  well  as  the  processes  in- 
volving the  combination  of  various  distributions  were  so  numerous  and 
complex  that  the  use  of  a  large  electronic  computer  was  mandatory. 

The  calculation  of  transmission  loss  values,     L     ,    was  accom- 

m 

plished  by  modifying  existing  programs  to  include  the  effect  of  the 

transmitting  antenna  only  as  a  function  of  distance  and  aircraft  height. 

The  resulting  propagation  curves  were  reduced  to  a  set  of  tables  for 

each  aircraft  altitude  that  included    L        information  for  each  of  the 

m 

three  propagation  models  discussed  in  Appendix  A.     Additional  tables 
for  distribution  functions  of  short-term  variability  and  aircraft  antenna 
power  gain  ratios  were  developed.     This  information  was  then  used  in 
the  Interference  Prediction  Program  to  generate  tables  of  interference 
ratios  for  various  values  of  time  availability  as  a  function  of  interfer- 
ence type,    station  separation  and  aircraft  location.     Values  of  normal- 
ized   D/U(95)    plotted  in  Section  5  were  extracted  from  these  tables. 

B.  1    Interference  Prediction  Program 

The  Interference  Prediction  Program  consisted  primarily  of  a 
series  of  loops  by  which  the  calculations  are  sequenced  and  systematic 
changes  of  parameters  achieved.     Figure  12  is  a  flow  diagram  of  the 
program  showing  the  series  of  loops  and  the  computational  sequence. 

The  primary  loop  considers  different  aircraft  altitudes.     At  the 
start  of  each  pass  around  the  loop  a  set  of    L        tables  (as  explained 
above)  for  a  particular  height  is  read  in.     This  loop  continues  until  the 
supply  of  tables  is  exhausted.     Altitudes  of  1,000,    6,  250,    12,000,    and 
18, 000  feet  were  used. 
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Figure   12,  Flow  Diagram  for  Interference  Prediction  Program 
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A  secondary  loop  considers  the  distance  along  the  great  circle 
path  on  the  ground  from  the  unde sired  station  to  the  aircraft,     d    .     An 
initial  value  of  5  nautical  miles  is  used  and  this  is  incremented  by  5 
miles  until  325  nautical  miles  is  reached. 

For  a  given  unde  sired  distance  an  inner  loop  is  used  to  calculate 
the  time  distributions  of  transmission  loss  for  the  two  unde  sired 
systems.     The  expression  for  the  cumulative  distribution  of  transmission 
loss  as  previously  given  is 

L(p.d)     =    Lm  -  [V(p,d) i*  VF(p,9)],  (27) 

Here  the  term    V(p,  d)    is  the  cumulative  distribution  representing 
long-term  variability  of  hourly  medians.     V    (p,  9)    represents  the  short- 
term  (within-the -hour)  variability. 

The  symbol    *    in  (27)  denotes  the  statistical  convolution  of  cumu-* 
lative  distributions.     Basically,    if  two  variables    x    and    y    are  statisti- 
cally independent  and  given  by  their  cumulative  distributions,    their 
statistical  convolution  is  the  cumulative  distribution  of  either  the  vari- 
able   z  =  x  +  y,    or  the  variable    z'   =  x  -  y,   depending  on  the  application 
of  the  result,    and  determined  in  the  machine  program  by  an  appropriate 
input  code.     The  cumulative  distribution  of    z    is  obtained  most  simply 
by  selecting    n    equally  spaced  percentage  values  from  the  individual 

distributions  of  both    x    and    y,    calculating  all  possible  sums     z,    =  x.  +  y., 

k         l         j 

and  forming  the  cumulative  distribution  of  all  values  of    z       obtained  in 

this  manner.     Similarly,    the  cumulative  distribution  of    z'     is  obtained 

by  forming  all  possible  differences     z    '    =  x.    -  y.. 

k  i       7j 

Successive  convolutions  are  handled  by  first  convoluting  two  dis- 
tributions and  then  using  this  intermediate  result  in  a  convolution  with 
the  third  individual  distribution.     The  convolution  process  is  cumulative; 
i.e.  ,    the  order  in  which  the  individual  distributions  are  used  does  not 
matter. 
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The  distributions  for  the  two  unde sired  systems  are  stored  and  a 

new  secondary  loop  which  considers  the  station  separation,     S,    is  now 

entered.     The  initial  separation  is  30  nautical  miles  and  is  incremented 

by  20  nautical  miles  until  a  separation  of  330  nautical  miles  is  reached. 

For  a  given  station  separation  the  desired  distance    d       is  now  defined 

as    d_  =  S  -  dTT    and  the  time  distribution  of  transmission  loss  for  the 
D  U 

desired  ILS  is  determined  in  the  same  manner  as  the  undesired  systems 
using  the  median  transmission  loss  value  from  the  appropriate  table. 
The  resultant  distribution,     L(p,d    ),   is  then  used  in  conjunction 
with  the    L(p,d    )    distributions  previously  stored  to  calculate    N{D/U(p)} 
for  the  two  desired -unde  sired  station  combinations  by  using  the  following 
equation: 

N{D/U(p)}    =    L(p,du)  *  [-  L(p,dD)]   *RA(p)  (28) 

where    R.(p),    as  explained  earlier,    represents  the  variability  introduced 
by  the  aircraft  antenna.     From  these  final  distributions  a  number  of 
percentage  values  between    p  =   1%    and    p  =  99%      are  selected  and 
punched  out. 

The  calculations  are  primarily  handled  by  three  subroutines: 
FALTG,   VPD,   and  INTERP.     This  technique  enables  the  program  to  be 
more  easily  changed  for  other  combinations  of  systems  not  being  con- 
sidered in  this  study. 

The  convolution  of  distributions  is  handled  by  the  subroutine 
FALTG. 

The  subroutine  VPD  calculates  the  predicted  long-term  variability 
(or  power  fading)  as  a  function  of  the  effective  distance  using  the 
methods  developed  for  a  continental  climate  and  presented  by  Rice,    et 
al.    [1965]  . 

With  the  extensive  use  of  tables,    the  need  to  interpolate  arises 
often,    and  linear  interpolation  will  not  always  suffice.     Therefore,    a 
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subroutine  INTERP  was  included  in  the  program  which  calculates  the 
n       Lagrangian  approximation  [Lance,    I960]  ,    where     n  <   m  -  1    with 
m    being  the  number  of  entries  in  the  table.     The  special  case     n  =  1 
results  in  linear  interpolation.     Asa  practical  matter,     n    is  usually- 
required  to  be  less  than  five.     This  is  necessitated  by  the  fact  that 
Lagrangian  interpolation  may  exhibit  large  oscillations  between  tabular 
entries  whenever  they  do  not  lie  on  rather  smooth  analytic  functions. 

It  should  be  realized  that  this  discussion  is  abbreviated,    and  that 
the  calculation  procedure  was  more  involved  than  indicated.     In  parti- 
cular,   the  method  used  to  extend  calculations  to  a  station  separation  of 
20  nautical  miles  was  not  discussed. 
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